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ABSTRACT 

We measure the shape of the velocity ellipsoid in two late-type spiral galaxies (Hubble 
types Sc and Scd) and combine these results with our previous analyses of six early- 
type spirals (Sa to Sbc) to probe the relation between galaxy morphology and the 
ratio of the vertical and radial dispersions. We confirm at much higher significance 
(99.9 percent) our prior detection of a tight correlation between these quantities. We 
explore the trends of the magnitude and shape of the velocity ellipsoid axes with galaxy 
properties (colour, gas surface mass density, and spiral arm structure). The observed 
relationships allow for an observational identification of the radial and vertical disk 
heating agents in external disk galaxies. 

Key words: galaxies: kinematics and dynamics - galaxies: individual: NGC 2280, 
NGC 3810 - galaxies: structure - galaxies: fundamental parameters 



1 INTRODUCTION 

More than half of the stellar mass in the local Universe 
is observed to resid e in disk galaxies l|Driver et al.l 120071 : 
IWeinzirl et al.| [2009'). yet the evolutionary history of these 
systems remains poorly understood. Simulations of galaxy 
formation tend to produce spiral galaxies with most of their 
stellar mass fraction in a bulge component, and a direct com- 
parison between observa tions and these simulati ons is usu- 
ally not straightforward (|Scannapieco et al.ll201QT ) . However, 
increasingly powerful simulations of disk galaxy formation 
are reaching a stage where internal dynamical properties, 
such as rotational and random velociti es, can also be in- 
vesti gated and compared to observations (|Scannapieco et all 
l201lh . These velocity distribution functions are fossil relics 
of the history of gravitational forces acting on stars and 
therefore of the evolutionary histories of spiral galaxies; for 
this reason, measurement of the stellar velocity distributions 
in local spiral galaxies is an important tool for probing the 
evolution of these systems. 

From observations in the solar neighbourhood it has 
been shown that the random motions of thin di sk stars 
corre l ates with their a ges, the age-er relation l|Wielenl 
1 19771 : iBinnev et alj|200d l. This relation is most commonly 
parametrized as a power law of the form u a T, with ex- 
isting measurements of a span ning the range of 0.2 — 0.5 
(see e.g. iNordstrom et al1l2004l for results from the Geneva- 



* Based on observations made with ESO Telescopes at the La 
Silla Observatory under programmes 074.B-0550(A) & 078. B- 
0152(A). 

f E-mail: jgerssen@aip.de 



Copenhagen survey). However, the age-a relation has also 
been interpreted as saturating at a constant dispersion af- 
ter ~ 5 Gyr (e.g. ICarlberg et aTll 19851: ISeabroke fc Gilmord 
120071 : ISoubiran et al.l [20081) or as several discr ete increases 
in velocity dispersion ( Quillen fc Garne"ttll200ll ) . As a result 
of the observational uncertainty, it is not straightforward to 
identify the disk heating mechanism responsible for the rise 
in velocity dispersions. 

Additional constraints on the disk heating mechanism 
are provided by the three-dimensional distribution of stellar 
velocity dispersions, whose magnitude and direction desig- 
nate the axes of an ellipsoid; th is trivariate Gaussian func- 
tion was originally proposed by ISchwarzschiidl (|l907h - The 
so-called velocity ellipsoid is parameterized by its two axis 
ratios, o^/ctr and ct z /o\r. In an axisymmetric disk with stel- 
lar orbits not too far from circular (the epicycle approxima- 
tion), o^/o-r depends only on the circular velocity and not 
on any disk heating mechanism. Measurements of the ratio 
o~ z /<7r can thus be used to probe the velocity distribution 
and therefore to constrain the heating processes in galactic 
disks. 

In the Milky Way, the age-cr relation and veloc- 
ity ellipsoid shape (u z /an = 0.5) have been stud- 
ied extensively, and a number of heating mechanisms 
have been proposed to simultaneously explain the ob- 
servations. These in clude encounters with giant molecu- 
lar c l ouds (GMCs. ISpitzer fc Schwarzschiid Il95ll : iLacevi 
1 19841 : lLacev fc Ostrikerl 1 19851 )." p erturbations from irregu- 
lar and transient spiral s t ructurelBarbanis fc Woltierl 1967 ; 
Jenkins fc Binnevl Il99p| : iFuchs! (200 it iMinchev fc Quillenl 



20061) . perturbations from stellar bars ( Saha et al 
dissolution of young stellar clusters (jKroupa 



2010), 



2002), 
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scattering by dark halo objects or globular c lusters 
l|Hanninen fe Flvnnl |2002| ; IVande Putte et al.l |200Sh . and 
disturbances by satellit e galaxies or minor mergers 
I Velazquez fc White! Il999h . This plethora of disk heating 
theories reflects the limited constraints provided by mea- 
surements in a single galaxy. However, many of these theo- 
ries are implicitly dependent on galaxy morphology; progress 
towards an understanding of secular evolution in disk galax- 
ies therefore critically depends on data from external galax- 
ies spanning a range of Hu bble types. 

The numerical study of ljenkins fc Binnevl jl990) is one 
of the few that can be directly compared with observations 
of disk heating in external galaxies. This theory combines 
two disk heating mechanisms, spiral transients and GMCs, 
in what appears to be the simplest explanation that is con- 
sistent with most of the observations. By varying the rel- 
ative importance of the two mechanisms they showed that 
disk heating, as quantified by the ratio a z /aR, decreases 
with the increase of s piral strength, e.g. along the Hubble 
sequence (but see also lSellwoodl l2008h . 

We have developed a technique l|Gerssen et al.l Il997f ) 
to measure the velocity dispersion ratio in galaxies from 
spectra obtained along the major and minor axis. We ap- 
plied this technique to six disk galaxies spanning a range 
in Hubble type from Sa to Sbc (jGerssen et all 1 19971 . 1200(1 
IShapiro et al.| [2003). In each system we were able to con- 
strain <j z /<jr by averaging measurements over radii in the 
range of one to two disk scale lengths. While the uncertain- 
ties on these ratios are large for individual galaxies, there 
appears to be a trend with Hubble type t hat is consistent 
with the predictions of ljenkins fc Binnevl |l990). However, 
the appl icability of this trend to later-type systems remained 
unclear (|Shapiro et al.ll2003l ). 

In this paper, we address this issue by directly mea- 
suring the velocity ellipsoid ratios in two late-type disk 
galaxies: NGC 3810 (Sc) and NGC 2280 (Scd). Section [2] 
presents our observations of these galaxies. Section [3] de- 
scribes our kinematic extraction and modeling procedure. 
Section [4] presents our derived velocity ellipsoid ratios and 
combines these results with our previous work on early-type 
disks to examine the ratio of vertical to radial velocity dis- 
persions across the Hubble sequence. The new results are 
consistent with the trend observed earlier and suggest that 
the relative effectiveness of vertical to radial secular evolu- 
tion processes diminishes as a function of galaxy morphol- 
ogy. We additionally combine our data with literature mea- 
surement s of the velocity ellipso id in elliptical and lenticular 
galaxies (|Cappellari et ahll2007h to reveal intriguing trends 
across all galaxy morphologies. In Section [5j we use our re- 
sults to identify heating agents at work in spiral galaxies, 
and in Section [6] we conclude. 



2 OBSERVATIONS & DATA REDUCTION 
2.1 Spectroscopy 

To probe the shape of the velocity ellipsoid at the end of 
the Hubble sequence, we obtained long-slit spectra along the 
major and minor axes of two late-type spiral galaxies. The 
galaxy properties are listed in Table [T] and the observing log 
is given in Table[2] All spectra were obtained with the EMMI 



Table 1. Galaxy properties 



Name 


Type 


Rcdshift 


i 


Scale length 


Reference 






(km s -1 ) 


(deg) 


(arcsec) 




NGC 2280 


Scd 


1899 


64 


19 


1, 2 


NGC 3810 


Sc 


993 


32 


16 


1, 3 



1) Redshifts and morphological classifications are from NED. 

2) Scale length and inclination are deri ved in section 12.21 

3) Scale length and inclination are from lKnapen et a! 1 120031) . 



Table 2. Observing log 



Galaxy 


Position angle 


Tea;p 


Date 




(deg) 


(min) 




NGC 2280 








Major axis 


-22 


156 


Feb 2005 


Minor axis 


-112 


160 


Nov 2006 


NGC 3810 








Major axis 


25 


210 


Feb 2005 


Minor axis 


115 


180 


Feb 2005 



spectrograph on the NTT during two runs in February 2005 
and November 2006. We used the red medium dispersion 
grating (#6) to study the wavelength range 4800 - 5500 
A, which covers key interstellar emission ([OIII], H/3) and 
stellar absorption (Mg b, several Fe) features. The slit length 
was 300", and the slit width was 1", resulting in a spectral 
resolution of 0.4 A (~ 23 km s _1 ). Measurements of the 
widths of unresolved arclines in our calibration spectra are 
consistent with this value. 

Observations were taken in exposures of either 1200s or 
2400s, which were interspersed with calibration exposures 
of a He+Ar lamp. In addition, we obtained spectra of sev- 
eral stars of spectral types GO III through K3 III, for use 
in extracting the stellar kinematics in our sample galaxies 
(i; 13. 1[) ; these stars are the dominant component in the Mgb 
and Fe stellar absorption features of spiral galaxies. To dis- 
tribute the light from the template stars relatively uniformly 
on the detector and to prevent overexposure, the telescope 
was defocused during the observations of these stars. 

Preliminary reduction steps (bias subtraction, flat- 
fielding) were done with our own IDL scripts, customized 
for EMMI data, and subsequent reduction steps (wavelength 
calibration, coadding individual exposures, extraction of 1- 
dimensional stellar spectra) were performed using standard 
IRAhQ packages. For the kinematic analysis, we interpolate 
the spectral axis to logarithmic units in wavelength (linear 
in velocity) and spatially bin each galaxy spectrum to a min- 
imum signal-to-noise (S/N) ratio of 10. 



1 IRAF is distributed by the National Optical Astronomy Ob- 
servatory, which is operated by the Association of Universities 
for Research in Astronomy (AURA) under cooperative agreement 
with the National Science Foundation. 
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2.2 Photometry 

Since the quantities of interest are properties of the galaxy 
disks and not of their bulges, we use broad-band images to 
measure the radii at which the disk dominates the stellar sur- 
face brightness. At these radii, the derived kinematics probe 
the disk stars and not the bulge and/or bar; we therefore 
limit our analysis (see § 13.21) t o thes e radii . 

For NGC 3810, iKnapen et ail |2003r i have measured 
the disk scale length to be 16" and the disk to dominate 
the galaxy light at r > 5". For NGC 22 80, we use an 
archival //-band image (jEskridge et al.ll2002r i. obtained with 
the CTIO 1.5m telescope and available via the NED extra- 
galactic database. We fit this image with the IRAF task 
ELLIPSE and derive an inclination of 64 degrees (from the 
best-fit axis ratios) and the surface brightness profile shown 
in Figure [T] The best-fit exponential disk model to this pro- 
file has a scale length of 19" and do minates the galaxy light 
at r > 10". (Using a .ft'-band image. iKassin et al. (|2006l ) de- 
rive a scale length of 27 arcsec.) The radial surface bright- 
ness profile (Figure 1) hints at the presence of a small bar 
in NGC 2280. Between the linear behaviour of the surface 
brightness profile at large radii (the exponential disk) and 
the inner bulge profile, there is a distinct third component. 
This component dominates the surface brightness profiles at 
radii of 5-10"; the combined shape of the surface brightness 
profile and fiat minor axis velocity dispersions (see § I4,1[I at 
these radii are consistent with the pr esence of a bar (e.g. 
iPrieto et al.1l200ll ; iGerssen et ai]|2003r ). 



3 ANALYSIS 

3.1 Stellar and Gas Kinematics 

For each binned spectrum, we extra ct the stellar kinematics 
using the ppxf method developed bv lCappellari fc Emselleml 
(|2004r ). We modeled the absorption line kinematics from our 
galaxy spectra as convolutions of a template stellar spec- 
trum with line-of-sight velocity distributions (LOSVD); the 
observation of the template stars with the same instrumen- 
tal set-up as the galaxies allows the instrumental broaden- 
ing to drop out during this fitting. Given the relatively low 
S/N of our data, we did not attempt to fit higher order 
moments of the LOSVD (e.g. h^,hi, etc) and instead as- 
sumed the LOSVD to be Gaussian. Since the outer disks of 
late-type galaxies are faint and the spectr a are dominated 
by sky emission, we employ the method of IWeiimans et all 
(2009), in which a median sky spectrum taken from the edge 
of our long-slit is included as a template in the fitting pro- 
cedure. We extract the stellar kinematics using several GUI 
and Kill templates. An optimal combination of these tem- 
plates did not significantly improve the results, so we chose 
the K1III giant HD114971, which yields marginally better 
X 2 , to extract the stellar kinematics in both galaxies. 

The derived stellar velocities and velocity dispersions 
along the major and minor axes of our sample galaxies are 
presented in Figure O In these plots, the abscissa is the 
radius in the plane of the disk, which for the minor axes 
was obtained through deprojection. The velocity dispersion 
profiles for both galaxies s how central drops. Th is is a fairly 
common phenomenon fe.g. lComeron et al.ll2008T l and is usu- 
ally attributed to a central disk or nuclear bar. These central 




40 60 
R (arcsec) 

Figure 1. The radial surface brightness profile of NGC 2280 
derived from an archival //-band imag e obtained with the CTIO 
1.5m telescope jEskridge et alj 120020 . The best-fit exponential 
model to the outer disk (solid line) has an exponential scale length 
of 19". The dashed line is a guide (not a fit) delineating the bulge 
component. The presence of a third component, perhaps a small 
bar, can be inferred from the excess of light between 5 and 10 
arcsec. 



velocity dispersion drops do not affect the results presented 
in this paper as their influence is confined to radii that we 
exclude from our analysis. 

To probe the circular velocity in our target galaxies we 
use the [OIII] emission line at 5007 A. We fit this line in 
each spatial bin with a Gaussian profile to derive the best- 
fit centroids. As the lines are unresolved in our data (width 
< 23 km s~ ) we do not attempt to quantify the gas velocity 
dispersions. In the following, we assume the measured gas 
velocities to be a good approximation of the circular velocity 
in our galaxies. 



3.2 Modeling 

To measure the velocity ellipsoid in our late-type galaxies, 
we model th e stellar and gas kin ematics using the method 
developed bv lGerssen et alj (|l997h and subse quently applied 
to a l arger range of morpholog ical types by IGerssen et al.l 
[|2000ri and lShapiro et all l|2003h . 

Briefly, the line-of-sight (los) kinematics along the ma- 
jor and minor axes of intermediate inclination galaxies probe 
different projections of the velocity ellipsoid, namely o§ and 
<jr along the major axis and (tr and ct 2 along the minor 
axis: 



"major = SU1 2 i + o\ COS 2 % 
2 2 • 2 • 2 2 • 

fj minor = (t q sin i + <j z cos t 



(1) 

(2) 



The epicycle approximation (e.g. Eq. 2 in lShapiro et al.l 

2003) for stellar orbits in a disk yields o^/an as a function 
of only the circular velocity V C {R) of the system. Assum- 
ing axisymmetry, measurements of <ti os along the major and 
minor axes and of V c along the major axis fully constrain 
the shape of the velocity ellipsoid. Including the stellar ro- 
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Figure 2. Observed kinematics of the two late-type spiral galaxies, left: NGC 2280 (type Scd) and right: NGC 3810 (type Sc). In all 
panels the kinematic information is given in km s" 1 . Filled circles represent major and minor axis stellar data; open circles represent 
gas data. Error bars for the stellar velocities are shown but are generally smaller than the plot symbols. The top panel for each galaxy 
shows the gas major axis rotation curve (blue open circles), the stellar major axis rotation curve (black filled circles), and the stellar 
minor axis velocities overplotted (red filled circles, roughly constant at km s — 1 ). The galaxy systemic velocities have been subtracted. 
The lower two panels for each galaxy show the major and minor axis stellar velocity dispersions respectively. The abscissa is the radius 
in the disk, which, for the minor axis, was obtained through deprojection. 



tation curve V(R) and modeling the asymmetric drift as a 
function of only the disk scale length a nd or (i.e. assuming 



Table 3. Best-fit model parameters 



no tilt orz to the velocity ellipsoid; see iGerssen et al. 1 119971 ) 
then overdetermines the system and provides an additional 
constraint to the model: 



v 2 = V? 



2 



R d 2 1 R 8V C 
--R— HffR )-- + — — 



(3) 



For a more det ailed descript i on of this method the 
reader is referred to IShapiro et all (|2003l ). In practice, the 
low intrinsic luminosity and small velocity dispersion of our 
lat e-type galaxies (compared t o the early-type disk s studied 
bv lGerssen et al.lll997l . l200d : IShapiro et afll2003l ) resulted 
in only a few data points sampling the radial behavior of 
the kinematic profiles. This paucity of data points forced us 
to include the emission line kinematics directly in the mod- 
elling rather than use them as an independent check on the 
best-fit results in our previous analyses. 

The four observables (cr ma j or , a m i nor , V c , and V) are si- 
multaneously fit to produce a z (R), ur(R), and V C (R). The 
velocity dispersions as a function of radius are assumed to 
be exponentially declining with identical scalelength hkin, 
parametrized by cr Zt o, (Jr^, and hki n , and the circular ve- 
locity is assumed to be a power law with radius, described 
by V c fi and power-law index a, for a total of five fitted pa- 
rameters. 

As our data set is rather small with fairly large error 
bars on each point, we employ two parameter space mini- 
mization techniques: the Levenberg-Marquardt method and 
a simulated annealing method. In both cases we follow the 



Parameter 


NGC 2280 


NGC 3810 


y c , (km s- 1 ) 


92 ±21 


142 ± 15 


a 


0.26 ±0.07 


0.11 ± 0.04 


(Tfl,o (km s _1 ) 


126 ± 29 


183 ± 24 


c z ,o (km s" 1 ) 


31 ±24 


54 ±22 


ft-ltin (arcsec) 


54 ±48 


22 ±3 


CTz/cTR 


0.25 ±0.20 


0.29 ± 0.12 



implem entation of these method as described in lPress et al.l 
(| 19921 ) . Both methods fit our data equally well in a \ 2 sense. 
The resulting best-fit models are overplotted on the data in 
Figure [3] and described in Table [3] 

The complexity of the dynamical modeling renders di- 
rect error propagation not straightforward in our analysis. 
We estimate the confidence limits using 500 Monte-Carlo 
realizations of the original data set. We analyze each real- 
ization using both \ 2 minimization techniques and find that 
they yield similar estimates of the uncertainties. 



4 RESULTS 
4.1 NGC 2280 

The Scd galaxy NGC 2280 is the latest Hubble type in 
our full sample of galaxies with measured velocity ellipsoid 
axis ratios. Kinematic data at radii less than 10 arcsec are 
dominated by the bulge and bar components in this system 
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Figure 3. Similar to Fig. 2 but showing the observed kinematic quantities at absolute radii and excluding data contaminated by bulge 
light, left: NGC 2280 and right: NGC 3810. All kinematic quantities (filled circles: stars, open circles: gas) are in km s — 1 . The best-fit 
model s to the observed ste llar kinematics (Eq. 3, 1 and 2 respectively from top to bottom) are shown as solid lines. Unlike our previous 
study jShapiro et alj |2003). where the gas velocities were used as a consistency check on the circular velocities (dotted lines), we include 
the gas velocities directly in our fitting procedure to better constrain the velocity ellipsoids in these later- type disk galaxies. 



and are therefore excluded from the analysis. The remain- 
ing dataset is limited in its extent but sufficient to derive 
a best-fit model, albeit with large uncertainties. The kine- 
matic scale length parameter, in particular, has a large un- 
certainty. However, this parameter drops out of the primary 
quantity of interest, the ratio of the vertical to radial ve- 
locity dispersions. The best-fit model velocity ellipsoid axis 
ratio is a z /a R = 0.25 ± 0.20. 

4.2 NGC 3810 

The photometric data for NGC 3810 (see H2.2p indicates that 
the disk light should dominate beyond R > 5 arcsec in the 
Sc galaxy NGC 3810. However, at these radii the measured 
velocity dispersions are not consistent with a stellar disk but 
instead hint at the presence of a (hidden) stellar bar. This is 
particularly apparent along the major axis at negative radii 
where the velocity dispersion profile is flat out to 8 arcsec. 
To avoid the contaminating effect of a potential stellar bar 
on velocity ellipsoid modelling we conservatively exclude all 
kinematic data at radii less than 10 arcsec. While this limits 
the dataset, its final size and quality are similar to that of 
the NGC 2280 dataset. The best-fit model velocity ellipsoid 
axis ratio is <7 z /<jr = 0.29 ± 0.12. 

4.3 The Velocity Ellipsoid in Spiral Galaxies 

Combining these measurements of <t z /(tr for late- 
type galaxies with our r esults for earlier-t ype spirals 
(Cerssen et all Il997l , l200d : IShapiro et ail |2003h . we probe 
the shape of the velocity ellipsoid over the full extent of the 



Hubble sequence of disk galaxies. The results are shown in 
Figure [4] and reveal a clear trend in <t z /(Tr with morphologi- 
cal type. The combined data are best fit by a line with slope 
—0.12 ±0.05, with a 0.1% probability of no correlation. The 
addition of our new results in late-type galaxies confirm and 
strengthen the trend iden tified at much lower s ignificance by 
iGerssen et ail (|2000l ) and lShapiro et all (|2003h . 

We include in the data and linear regression in Figure [4] 
the Hipparcos value for th e velocity ellipsoid in th e solar 
neighbourhood (0.53 ± 0.07. iDehnen fe Binnevlll998l ) . which 
is consistent with results b ased on the RAVE data (0.55 , 
ICasetti-Dinescu et alj|201ll . but see also lVeltz et all 120081 ). 
This measurement is obtained at the solar radius (at ~ 2 — 3 
scale lengths in the Milky Way) and is therefore comparable 
to our measurements in external galaxies, radially averaged 
in each galaxy over a typical range of ~ 1 — 3 disk scale 

lengths. 

Using a different technique, Ivan der Kruit k, de Griisl 

1 19991 ) did not observe a clear trend in velocity ellipsoid 
ratio in galaxies of Hubble type Sb to Sd. However, their re- 
sults are statistical (the rms of determinations for multiple 
edge-on galaxies of the same type) and based on photome- 
try only. An advance of this technique combines photome- 
try and spectra of edge-on galaxies to measure the quantity 
sj Mj L(aR/a z )~ 1 ; assumptions of universal mass-to- light 
(M/L) ratios and/or Toomre Q parameters in late- type spi- 
rals can then be us ed to estimate the velocity ellipsoid ratio 
jKreeel et al.ll2005l ). For their sample of 15 edge-on galaxies 
of Hubble types Sb to Scd, these authors obtained a range of 
estimates of cf z //jr consistent with our direct measurements 
and following a marginal trend with Hubble type. 

A promising step forward is the use of integral field 
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Figure 4. Velocity ellipsoid ratio a z /(jR as a function of galactic 
type (Hubble stage T), The solid points s how the results that 
we obtained previously iShapiro et alj2003l) . The results that we 
derive in this paper for two late type spirals are shown as the 
open circles. Horizontal errors represent the uncertainty inherent 
in galaxy classification IjNaim et al,lll995T ). The filled square is the 
value in the Solar Neighbourhood derived from Hipparcos data 
toehnen fc Binne^ll99Sl) . 



Figure 5. Velocity ellipsoid ratio <j z /o-r as a function of 
inclination- and extinction-corrected galaxy colour (available in 
HyperLeda for 7 of our 8 spirals and the fast-rotating E/SOs). The 
black points ind icate our spiral g a laxy d ata, and the red squares 
are the data of ICappcllari et al.l 1120071 ) for fast-rotating E and 
SO galaxies. The linear fit (dashed line) is to the spiral galaxies 
only. The probability of no correlation, null hypothesis Ho, is one 
percent. 



units (IFUs) to observe velocity dispersions in disk galaxies. 
This is observationally more efficient than obtaining long-slit 
spectra, one at a time, along two (or more) position angles. 
Additionally, since IFUs uniformly sample velocity disper- 
sions along both azimuth and radius across the disk, the 
assumption of the epicycle appr oximation, employed here 
for long-slit data, can be relaxed. iNoordermeer et al.l (2008) 
use the PPAK IFU to explore velocity dispersions in disk 
galaxies and to constrain the shape of the velocity ellipsoid 
in one system, NGC 2985. Their measurement of o z /or 
~ 0.7 is consist ent with our previo us result for this sys- 
tem (0.75 ± 0.09; iGerssen et alj|2000t ) and confirms that the 
epicycle theory we assume is indeed applicabl e in this galaxy. 

In a s eries of conference proceedings IWestfall et al.l 
(2008, 2010) present their ongoing work with the SparsePAK 
and PPAK IFUs to constrain the velocity ellipsoid shape 
and its radial variation in a sample of spiral galaxies. These 
authors find a strong dependence on their modelled veloc- 
ity ellipsoid shape with their measurement techniques and 
assumptions; they present and describe an analysis compa- 
rable to ours for a single galaxy, NGC 3982 (Sb), for which 
they find <j z /or = 0.31 — 0.73 over a radial range of 1 — 2 
photometric scalelengths, broadly consistent with our results 
for galaxies of similar Hubble type. These authors also find 
some evidence for variation in the velocity ellipsoid ratio 
with radius, indicating the potential of IFU data in future 
studies of the velocity ellipsoid in external galaxies. 

4.4 The Velocity Ellipsoid across the Hubble 
Sequence 

The SAURON team has used the IFU of the same name 
to study the kinematics of elliptical and lenticular galaxies 
and have shown that many early-type galaxies, the so-called 



"fast rotators," are bulge-dominated gala xies that neverthe- 
less contain a significant disk component l|Kuntschner et all 
12009 ; iKrainovic et al.ll2008h . It is therefore interesting to 
investigate how the disks of these early-type galaxies are 
rela ted to those in sp ir al gala xies. 

ICappellari et all (|2007l ) have used axisymmetric 
Schwarzschild dynamical models to extract the three- 
dimensional orbital structure of a subsample of the 
SAURON galaxies and measure the shape of their velocity 
ellipsoids. These anisotropy measurements are luminosity- 
weighted, giving more weight to the high-density equatorial 
plane, and volume- averaged, giving more weight to larger 
radii; as a result, the global anisotropies (see table 2 of 
ICappellari et al.1l2007l ) are dominated by the disks and are 
therefore comparable to our measurements in later-type 
galaxies (Cappellari et al. private communication). Since 
Hubble T-type is less meaningful in early-type galaxies 
than in spirals, we cannot add these points directly to 
Figure 0] Instead we use galaxy B — V colour as a proxy for 
morphological type and plot this quantity against <t z /<jr in 
Figure [5] for both the SAURON sample and our sample of 
spiral galaxies. Colours are taken from HyperLedc0 and are 
listed for our sample in Table [3] 

The combined data span the Hubble sequence from E 
to Scd and show the strong correlation between velocity el- 
lipsoid ratio and galaxy colour, as expected from Figure 3J 
given the known relationship between galaxy colour and 
Hubble type. In Figure [5] we find a continuous trend of in- 
creasing anisotropy in bluer galaxies. Moreover, we find that 
the anisotropies of the E/S0 galaxies overlap with those of 
the earliest-type spirals in our sample. However, in the bluest 

2 http://leda.univ-lyonl.fr/ 
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range (B — V ~ 0.7 — 0.8) of the overlapping data, there is 
some evidence that elliptical and lenticular fast rotators are 
preferentially more isotropic than their spiral cousins. 



5 DISCUSSION 

5.1 Dependence of the Velocity Ellipsoid Shape 
on Galaxy Properties 

Figures [4] and [5] point to a deep connection between galaxy- 
wide properties and axial ratio of the stellar velocity ellipsoid 
in the galaxy disk. In this and the following sections, we 
examine the source of these relationships. 

The tighter of the two correlations is that between 
galaxy morphology, parameterized through Hubble T-type, 
and the shape of the velocity ellipsoid (Figure [4]). However, 
as a semi-quantitative descriptor of galaxies, T-type is a 
somewhat subjective property based on several qualitative 
criteria that are largely but not perfectly correlated. In the 
de Vaucouleurs RC3 classification used here, these criteria 
are the bulge-to-disk ratio and the tightness of the spiral 
structure. Moreover, a number of additional galaxy prop- 
erties not directly involved in the determination of T-type 
also correlate with morphology, such that earlier-type spi- 
rals tend to have higher masses, higher stellar mass sur- 
face densities, redder colours, and lower gas fractions than 
thei r late-type counterparts , albeit with large scatter (see 
e.g. iRoberts fc HavnesHl994T ). 

It is therefore not straightforward to determine which 
of the fundamental galaxy properties that correlate with T- 
type drive the relation in Figures [4] and [5] Broadly, these 
properties fall into two general categories: (1) structural 
properties related to the large-scale evolutionary history of 
a galaxy such as galaxy mass and relative prominence of 
the bulge, and (2) properties related to past, current, and 
future star formation processes such as galaxy colour, tight- 
ness of the spiral structure, and gas fraction. The interesting 
result from Figure [3] is that one or both of these aspects of 
a galaxy's evolution is intimately related to the heating of 
its disk. 

Moreover, Figure [5] reveals that these same processes 
are also at work in fast-rotating early-type galaxies. The 
more isotropic shape of the velocity ellipsoid in early-types, 
relative to late-types with similar colours, may indicate a 
difference in disk heating mechanisms at work in early- and 
late-type galaxies. One possibility is a difference in strength 
of a particular mechanism; another is that an additional 
heating agent may need to be invoked. 

5.2 Disentangling Disk Heating Agents 

Of the numerous theoretically studied disk heating agents 
(see § [TJ , each is expected to behave differently in the verti- 
cal and radial directions as a function of galaxy morphology. 
To constrain the nature of these agents in our spiral galax- 
ies, we therefore examine the individual axes of the velocity 
ellipsoid in our galaxy sample. (The evolution of the velocity 
ellipsoid in the SAURON sample of elliptical and lenticular 
galaxies is an equally interesting problem beyond the scope 
of this paper, and we do not discuss it further.) 

Using the exponential parametrization of a z and ctr 



with radius (§ 13. 2[) . we evaluate the magnitude of the 
dispersion at half the measured kinematic scale length, 
hki n /2. This radius falls within the range constrained by 
the observations for all but one of our eight galaxies. 
The local isothermal appro ximation of stellar disks predicts 
that h kin /2 fts h v hot (e.g. Ivan der Kruit fc Freeman! 1 19861 : 
iMartinssonl 1201 IT ). Our measurements of hkin are broadly 
consistent with this prediction, albeit with large scatter and 
noticeable outliers. 

The relationship of the magnitudes of the dispersions at 
hkin/2 to each other and to the shape of the velocity ellip- 
soid is shown in Figure [5] In this figure, the tightest trend, 
with a 1% probability of no correlation, is between the mag- 
nitude of the vertical velocity dispersions a z and the shape 
of the velocity ellipsoid (t z /or\ in stark contrast, there is no 
correlation between o z /ctr and the magnitude of the radial 
velocity dispersions or. Although there is evidence for a re- 
lationship between a z and <jr, a constant ratio (i.e. linear fit) 
is inconsistent with our direct measurements of variations in 
&z/<?r among galaxies and cannot be the best description 
of the data. In the left panel of Figure [5] we therefore plot 
the prediction for the relationship between a z and an based 
on the observed linear correlation between a z and ct z /(Jr 
in the middle panel. The projection of the (t z -<j z /or cor- 
relation into the (a z ,au) parameter space is shown as the 
dashed line. The uncertainty in this projection is shown by 
the shaded area, bounded by the la uncertainties on the lin- 
ear fit and by the line where (t z /ctr= 1. Within the scatter, 
the predicted relationship between a z and or matches the 
data. 

These results are a powerful new tool in probing the 
nature of disk heating agents at work across the Hubble 
sequence, and they demonstrate the feasibility of observa- 
tionally disentangling disk heating agents in spiral galaxies. 
In particular, Figure [6] places two important constraints on 
disk heating mechanisms: 

(1) There is a disk heating agent that is responsible 
for increasing both vertical and radial dispersions, which we 
call the "three-dimensional agent." The trend of a z with 
or demonstrates that, in external galaxies, both dispersion 
components grow simultaneously although not at the same 
rate, as is also obser ved in the solar neighbourhood (e.g. 
iHolmberg et~aT1 12007I ). The tight correlation between a z 
and cr z /tJR additionally requires that the three-dimensional 
agent responsible for the magnitude of the vertical disper- 
sions also determine the shape of the velocity ellipsoid, such 
that larger vertical dispersions are coincident with more 
isotropic heating and smaller vertical dispersions are coin- 
cident with more radially anisotropic heating. The correla- 
tions in Figures [4] and [6] therefore imply that the strength 
of the three-dimensional agent in a galaxy should correlate 
strongly with a z , <j z /or, and the galaxy properties associ- 
ated with T-type and correlate weakly with or. 

(2) It is likely, given the observed scatter in the o z -ctr 
relation, that there is also a second agent involved in ra- 
dial heating, which we call the "radial agent." The lack of 
correlation between or and o z /<jr suggests that this heat- 
ing mechanism has no effect on the vertical dispersions (and 
therefore on the shape of the velocity ellipsoid) and operates 
only in the plane of the galaxy. Additionally, the lack of cor- 
relation between or and o"z/o\r implies that the presence or 
strength of the radial agent in a galaxy should correlate only 
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Figure 6. Correlations between the velocity ellipsoid magnitudes and shapes for all eight galaxies in our sample. Here and in the following 
figures, Ho denotes the probability of the null hypothesis, i.e. that there is no correlation. Left panel: The vertical ellipsoid magnitudes 
as a function of the radial component magnitude. The dashed line is the projection of the best-fit correlation in the middle panel, see 
text. The shaded area is bounded by the one-sigma uncertainty on the projected correlation and by the line cr z /an= 1. Middle panel: 
The shape of the velocity ellipsoid for each galaxy in our sample plotted as a function of the vertical dispersion magnitude. The results 
are tightly correlated, dashed line. Right panel: The velocity ellipsoid shape as a function of the radial dispersion magnitude. There 
appears to be no relation. 



Table 4. Properties of disk galaxies with velocity ellipsoid shape and magnitude measurements. 



Galaxy 


Hubble Type 


D (Mpc) 




(B-V)Tot 


S H2 (M /pc- 2 ) 


Arm Class 


NGC 2460 


Sa 


20 


0.83 ± 0.35 


0.91 




12 


NGC 2775 


Sa/Sab 


19 


<1.02 ± 0.11 


0.90 


33 ± 6 


3 


NGC 2985 


Sab 


18 


0.75 ± 0.09 


0.74 


21 ± 7 


3 


NGC 488 


Sb 


30 


0.70 ± 0.16 


0.87 


25 ± 6 


3 


NGC 1068 


Sb 


16 


0.58 ± 0.07 


0.74 


41 ± 14 


3 


NGC 4030 


Sbc 


21 


0.64 ± 0.28 




16 ± 3 


9 


NGC 3810 


Sc 


14 


0.29 ± 0.12 


0.58 


28 ± 6 


2 


NGC 2280 


Scd 


26 


0.25 ± 0.20 


0.60 




9 



1) Hubble classifications are from NED. 

2) Distances are derived from the redshift values listed in NED and assuming a Hubble constant of 73 km/s/Mpc. 

3) Colours are from HyperLeda. 

4) Arm Class values are from Elmcgrccn & Elmcgrccn (1987). 



with or, and not with a z , o z /(tr, or the galaxy properties 
associated with T-type. 

Two of the most commonly invoked disk heating agents 
are GMCs and transient spiral structure (§[1]), with the for- 
mer expected to operate in three dimensions and the latter 
only in the plane of the disk. In Sections 15.31 and 15.41 we in- 
vestigate whether these heating mechanisms are consistent 
with our observational constraints on the three-dimensional 
agent and radial agent, respectively. 



5.3 GMCs as a Heating Agent 

We estimate the molecular gas surface mass densities for six 
of the eight galaxies in our sample for which the CO models 
of lYoung et al.l (|l995l ') are available. We compute E_h- 2 in 
each system as the average surface density over the same 
range in radii used in our velocity ellipsoid measurements. 
The results, with rather large error bars, are tabulated in 
Table H 

In Figure [7] the velocity ellipsoid shape and magnitude 



are shown as a function of E_h 2 . If molecular gas is the dom- 
inant "three-dimensional agent", we expect the gas surface 
density to correlate with o z , or (weakly), o z /or, and T- 
type. The radial dispersion shown in the right panel are 
indeed consistent with the expected increase in dispersions 
with gas surface density; however, a trend is not evident in 
the vertical dispersions (middle panel) . The probability that 
neither radial nor vertical dispersions are correlated with 
surface density is 19% and 44%, respectively. Moreover, E# 2 
is also not correlated with the velocity ellipsoid shape ct z /(tr 
and is instead consistent with s cattering around a constant 
value. iRoberts fc Havneslll994l show that the surface den- 
sity of molecular gas increases with Hubble T-type, although 
this trend is not clearly reflected in our much smaller sample 
(not shown). Together, these relationships point to GMCs as 
only marginally consistent with being the three-dimensional 
agent. 

Theoretically, three-dimensional heating via GMCs has 
been shown to drive the velocity ellipsoid axis ratio from an 
initially isotropic distribution to a small equilibrium value 
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(e.g. I Jenkins fc Binnevlll990l ). Interestingly, ISellwoodl (|2008h 
argued that this mechanism alone can explain the observed 
disk heating in the Milky Way, with the degree of anisotropy 
in the distribution of the scattering events determining the 
final equilibrium ratio o z /or. However, our observed vari- 
ation in o z /or with T-type in external galaxies is larger 
than can be comfortably explained with GMCs as the only 
heating mechanism. 

The identification of GMCs as the three-dimensional 
heating agent is not straightforward. GMCs are consistent 
with the behavior required by the data b ut leave room 
for alternative explanations. Ijenkinsl (|l992T ) likewise notes 
that the observations of constant disk scale height with 
radius in external galaxies is difficult to reconcile with 
the observed exponential distribution of molecular clouds 
if GMCs are important disk heating agents. Other poten- 
tial three-dimensional agents include transient asymmet- 
ric spiral structure, dark halo objects (black holes and 
dark matter substruct ure), globu l ar clusters, and accret- 
ing s a tellite galaxies ijSaha et all 20ld; Laccv & O s triker! 



1985; Ha nninen fc Flvn 



Velazquez fc White! [199' 



'lvnn l 
19991 ) 



20021 ; IVande Putte et al.l 120091 ; 

Of these, numerical simula- 
tions have demonstrated that dark halo objects and glob- 
ular clusters are unlikely to exist in sufficient num- 
bers t o produce the obs e rved heating dLacev fe Ostrikeil 
19851; lHanninen fc Flynnl |2002| ; IVande Putte et al.l 1200a 



Benson et al.ll2004l ). leaving spiral structure and accreting 
satellites as viable mechanisms. It is not obvious how the 
accretion history of a galaxy might be quantified and sub- 
sequently compared to the observed shape of the velocity 
ellipsoid, so observational tests of the viability of this alter- 
native heating mechanism are not straightforward. In con- 
trast, spiral structure is a direct observable, and we examine 
the effects of this potential heating agent in the following 
section. 



5.4 Spiral Structure as a Heating Agent 

To investigate the role of spiral structure on disk heat- 
ing we use the arm class parameter, developed by 
lElmegreen fc Elmegreenl l| 19871 ). to quantify the strength 
and symmetry of the spiral structure in galaxies. They define 
this parameter to describe the orderliness of spiral structure 
from flocculent (class 1) to grand-design (class 12). The arm 
class parameters for our sample are listed in Table 0] 

We compare the velocity ellipsoid shape and magni- 
tudes in our sample galaxies to their arm classes in Fig- 
ure [5] and find a strong correlation between or and arm 
class. This can be interpreted as spiral structure preferen- 
tially for ming in Toomr e unstable disks (Toomre Q stability 
criterion; lToomrdll964h . The Toomre Q is proportional to 
or, such that larger velocity dispersions create more stable 
systems, which are then less res ponsive to the formation of 
(transient) spiral structure (e.g. lBinnevll2012l ). In our sam- 
ple, we indeed find the most clear spiral structure (class 12) 
in those galaxies with the lowest or. 

In Figure [5] there is a marginal relationship between 
o z and spiral structure at much lower significance than 
that with or, and no trend is evident between cj z /(jr and 
arm class. This latter result is consistent with the ab- 
sence of a trend betwee n arm class and T-type found by 
lElmegreen fc Elmegreenl l|l987l ) and confirmed in our much 



smaller sample (not shown) . Spiral structure thus correlates 
primarily with or, weakly (if at all) with o z , and not with 
o z /or nor T-type. Based on the requirements of different 
heating agents described in Section [5[2] spiral transients are 
excellent candidates for the "radial agent" . 



6 CONCLUSIONS 

In this paper, we present the kinematics and modelled ve- 
locity ellipsoid shape, i.e. the ratio of the vertical and radial 
velocity dispersions o z /or, in two late- type disk galaxies, 
NGC 2280 (Scd) and NGC 3810 (Sc). Combining these re- 
sults with our previous measurements of this ratio in six 
early type disks demonstrates that the shape of the velocity 
ellipsoid is strongly correlated with Hubble type and with 
galaxy color. Moreover, early-type galaxies (E/SO) with sig- 
nificant disk components appear to follow a related trend. 

Examination of the relationship of vertical and radial 
dispersions (o z and or) with each other and with o z /or 
shows that o z /or correlates strongly with o z but not at all 
with or. This suggests that kinematic observations alone 
can disentangle the effects of disk heating mechanisms that 
act in three-dimensions from the effects of those that act 
predominantly in the plane of the disk. 

Two commonly invoked disk heating mechanisms are 
stellar scattering off of giant molecular clouds (as the three- 
dimensional agent) and perturbations from transient spiral 
structure (as the radial agent). To test the applicability of 
these theories to our data across the full range of morpholog- 
ical types, we probe the relationship of molecular gas sur- 
face mass density and of spiral arm structure to our ob- 
served dispersions (o z and or) and their ratio (o z /or). The 
data show that spiral structure varies between galaxies ex- 
actly as expected of the radial heating agent, and there- 
fore it seems likely that spiral structure is responsible for 
(morphologically-independent) heating in the plane of galac- 
tic disks. The data are tentatively consistent with GMCs be- 
ing the three-dimensional agent, although the gap between 
the expected and observed variation between galaxies for 
the three-dimensional agent leaves room for alternative in- 
terpretations. 

These results suggest exciting possibilities for future 
work in directly probing disk heating agents in external 
galaxies. Open questions include the nature of the three- 
dimensional agent and its relationship to the molecular con- 
tent of galaxies, the possible relationship of stellar popula- 
tions to disk heating mechanisms, and the relationship of 
disk heating mechanisms in fast-rotating early-type galaxies 
to those of spirals. 
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Figure 7. The veloci ty ellipsoid shape a nd magnitudes as a function of the H2 gas surface density. Gas densities are estimated from the 
CO measurements of I Young et al,l il995l ) and are computed as average values over the radii used in our kinematic analysis. Left: The 
velocity ellipsoid shapes are not correlated with the molecular gas surface density. Middle: There is a hint that the vertical component 
of the velocity dispersion is correlated with 2^ 2 , but the scatter is too large to state this conclusively. Right: The radial component 
increases with molecular gas density. 
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Figure 8. The velocity ellipsoid shape and magnitudes as a function of arm class, as defined in lElmegreen &: Elmegreen! l|l987V ) to 
quantify the or derliness of spiral structure from fiocculent (class 1) to grand-design (class 12). Note that there is no arm-class 10 and 11 
(cf., Figure 1 in Elmcgrccn & Elmcgrecn 1987). Left: The velocity ellipsoid shape is not correlated with arm class. Middle: The vertical 
magnitude of the ellipsoid decreases with arm class. Right: There is a clear trend between arm class and radial component, as expected 
from the Toomre Q criterion, see text. Note that for plotting purposes we have added small offsets to galaxies with the same arm class. 
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